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ABSTRACT 

A theoretical model is developed for the wind velocity 
field in the atmospheric near-water layer above ocean waves, 

A turbulent exchange coefficient due to the wave-caused 
velocity fluctuations is introduced into the equations to 
account for the corresponding turbulent nature of the motion 
and the resulting energy fluxes. The equations of motion 
for air are derived in a manner to allow for variation of 
the turbulent exchange coefficient with height. A stream 
function is introduced and a numerical solution obtained by 
utilization of the Richtmyer method. The profiles of wave- 
caused Reynolds stresses and their resulting energy trans- 
fers are discussed and compared with observational data. 
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z - height 

2 
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I. INTRODUCTION 



Energy and momentum transfer to water waves are pri- 
marily dependent upon the conditions of the atmosphere 
immediately adjacent to the ocean surface. In the case 
of air flow above an open body of water the atmosphere- 
ocean interactions give rise to wave perturbations in the 
upper layer of the water and in the adjacent atmosphere. 

It has been suggested that the wave perturbation charac- 
teristics be used to determine the energy and momentum 
transfers across the air-sea boundary. 

Investigation of the nature of these energy and momen- 
tum transfers have been the purpose of numerous studies 
in recent years. Such studies of wind-wave coupling were 
initiated with the quas i — laminar model of Miles (1957 ) , 
who allowed the wave-induced motion to interact with the 
shear flow only at the critical level. Models progressed 
through the studies of Phillips (1957) and Benjamin (1958) 
to those more recently of Yefimov (1970) and Reynolds 
(1972), wherein there is also interaction with the turbu- 
lent Reynolds stresses. 

In the present study, wave related momentum transfers 
were investigated through the use of a numerical model. 
This study is similar to that of Yefimov (1970) , but it is 
more general and allows for dependence of the mean wind 
on time and variation of the turbulent coefficient with 
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height. This model utilizes an initial value numerical 
procedure similar to that of Newman (1969). 

In addition, specific data derived from the statisti- 
cal study of Davidson and Frank (1973) will be utilized 
to define solution parameters to generate the numerically 
derived stress patterns. 
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II. FORMULATION OF THE MODEL 



The basic of the problem is the equation of motion. 
Since the stress is expected to be independent of y, two- 
dimensional motion is assumed with the resulting equation 
in the form: 

4— + W • V W = - ~ VP + ^ , (2.1) 

0 L P 



where 



F x • (-S' 2 ) + and 



F z ■ 57 + li ■ 

Here x is the horizontal axis; z is the vertical axis; p 
is the density of air; P is the non-hydrostatic pressure; 
u f and w ! are the turbulent velocity fluctuations. The 
stationary part of each variable (time average) will be 
denoted with a straight line above it and the nonstation- 
ary part (space average) by a waved line. The stationary 
velocities do not include the random components which 
would disappear with a suitable time average. The veloc- 
ity (W) will be defined as: 



W = (u + U) U. + w Jk (2.2) 

/ J 

3 U 

where U = U(z) and the factor -r — will be neglected in the 

d t 

development of the equations. For simplicity the problem 
will be restricted to those cases of a single wave number. 
This neglects the effect of wave self-interaction which 
would affect other waves. 
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The turbulent velocity fluctuations, 



u f and w ! are 



assumed to have a dual nature. They are gener 
gradient of the mean velocity in the near wate 
also by the wave-caused perturbations of the v 
The tangential Reynolds stress (u’w 1 ) will be 
as : 



where u and 
of wave-cau 
the coeffic 
the generat 
ity gradien 
lent exchan 
fluctuation 



3U 



3u , 3w> 



T — —u 1 W f = K ~ + V 

3z 3z 3x 



w are the horizontal and vertical 
sed velocity perturbations. The fa 
ient of turbulent exchange which is 
ion of turbulent fluctuations by th 
t. The element V is the coefficien 
ge which is related to the wave-cau 
s. It will be assumed that V >> V 



m 

except in the thin laminar sublayer in which w 
gleet the effect of turbulent viscosity. Both 
cients, K and V, will be assumed to be constan 
x-direction and Vis a function of height. In 
the coefficient K will be assumed to be given 

K = K U* z 

which increases linearly with height and leads 
disturbed mean velocity profile 

u * z 

u ■ — in <— + • 



o 



ated by the 
r layer and 
elocity . 
expressed 

, ( 2 . 3 ) 

comp onen t s 
ctor K is 
related to 
e mean veloc- 
t of turbu- 
sed velocity 
everywhere 
e shall ne- 
the coeffi- 
t in the 
addition , 
by : 

to the un- 



( 2 . 4 ) 
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Incorporating the continuity equation with the Bous- 
sinesq approximation, (V*W) = 0, the equations may be 
written in component form as: 

9u .u/.n ^ 9U 1 9P . 9 e — 2> , 9 rir 9U . .3u , 3w. n , 0 

— + \V* Vu+w— — (-u T ) + x— [K~ 4* v (— + -r— ) ] , (2 . 5) 

dt 9z p 8x 9x 9z 9z 9z 9x 

9w .\u.T 7 1 9P 9 rir 9U a1 /9u 9w. , 9 , — 2 \ / 9 

7 T- + \V* Vw « - — — 4 - — [k— 4 * V (— 4* ) ] 4 * — (-w ) . (2.6) 

9z p 8z 9x 9z 8z 9x 9z 

The vorticity (£) is defined to be: 



9u 



9w 



^ 9 z 9x 



-2 ^ 3 —2 
2 



Normally u T “ = ■— w 1 “ but for this work the assumption of 



i 2 ,2 

u = w 



will be incorporated into the vorticity equation which may 
now be written in the form: 

3C . ... nr , 3W 3W , 3 , 31k 

77 + W ’ V?+ 37 ’ Vu ' 37 * Vw + ^ (w ^ } = 



3z v 3c' 



3 3 u 



3 z 



r V 3U . 3u 3Wv . r 

7 1k 3¥ +v <37 + * 



3 3 



w 



9 z 9 z 



9x 



] (2,7) 



If the continuity equation, 



+ 12 = o 

3x 3z 



is used, the vorticity equation becomes: 

J. IT a. ~ ^ 2l 

3 1 



9? , „ 3? , 3‘U , 3? , 

7T” + u - — 4* w 4* u ~ — 4* W 



9x 



9z ' 



9x 



17 = 

9z 



9 2 rtr9U , ,/9u , 9w xn / 9 3 u , 9 3 w 

— 7 IK 3l + v <37 + 17 + 77 

9z 9x9z9x 



) ( 2 . 8 ) 



The products of the perturbation quantities are neglected 
except where they affect the mean flow. This is consistent 
with the consideration of a single wave number. 
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The mean quantities are balanced separately and sub- 
tracted from (2.8). This results in the two-equation 



system : 



at 



+u 



3x 



+ W' 



3 2 U 



3 z ‘ 




tv ( 



3u 
3 z 



+ 



3w 

3x 



))-V[ 



3 3 u 
3x 2 3 z 



+ 



3 3 w j 
3x 3 



(2.9) 



u f + „ K. 2i_ 

3x 3z rx 2 

o z 



1=) 



( 2 . 10 ) 



With the introduction of stream function notation, the 
velocities can be defined as: 



which 




u 



t ransf onus 



+ \p ) = 
zz r xx 



3j l 

3 z 
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(2.9) and (2.10) into: 

-U ( Tb + rp ) + to U 

r XZZ T xxx r x zz 



( 2 . 11 ) 



+ - t v (i|> - \p )] - V [ip - \p ] , (2.12) 

~ 2 r zz xx r xxzz r xxxx 1 

3 z 

^ ~ 

\p \p + ip \p - \p \p - ip~\p = K\p (2.13], 

z r xxx r z r xzz T x zxx r x zzz r zzzz 

Yefimov (1970) simplified the solution by assuming K and 
to be of the same order. This assumption will be incor- 
porated into equation (2.13) by replacing K by V. 

The analytic solution for the special case, U = 0, to 
(2,12) and (2.13) has been considered by Yefimov and 
Pososhkov (1970). Their results have been adapted and in- 
corporated as a test of the numerical technique as outlined 
in Appendix C. 
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The stream function will be defined as: 
ip = ^ Q (z,t)+ A(z , t) cos [k(x-ct) ]+ B (z , t) sin[k(x-ct) ] (2.14) 

where A(z,t) and B(z,t) are the stream function amplitudes. 

Substitute the expressions for the various derivations 
of \p into (2.12) and (2.13) and equate coefficients of the 
cosine terms to produce: 



f-(A - Ak 2 ) = (U-c)[Bk 3 - B k] + BkU + 

9 1 zz zz zz 



9 z 



— [ V ( A + Ak )] + V [A k + Ak ] 

2 Z Z ZZ 



The coefficients of the sine terms are: 

- Bk 2 ) = -<U-c) [Ak 3 - A 22 k) - AkU 22 + 



zz 



— [ V ( B + Bk 2 ) ] + V [B k 2 + Bk^] 
2 zz zz 



(2.15) 



(2.16) 



The mean velocity equation is 



|-[A B -A B-AB+AB ] = (2 . 1 7 ) 

2 l z zz zzz zz z zzz ozzzz 

Equations (2.15) and (2.16) determine the wave-caused 
velocity perturbations and therefore are the basic equa- 
tions for solving the problem. Equation (2.17) is used 
for the calculation of the mean velocity perturbation u. 
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III. SOLUTION OF THE PROBLEM 



A. INITIALIZATION 

To initiate the solution of equations (2.15) and (2.16), 
the undisturbed mean velocity was assumed to be in the 
form: 

U * 

U ( z ) = — In (— + 1) , (3.1) 

o 

as shown in Figure 1. 

The functions A(z,t) and B(z,t) were initialized 
through : 

A(z,0) = A(0,0) sinh [k (Hl-z ) ] / s inh (kHl) , (3.2) 

and B ( z , 0 ) - 0 , (3.3) 

where H is the height of the boundary layer and HI the 
height one grid point less than H. 

B. BOUNDARY CONDITIONS 

Yefimov (1970) defined the boundary conditions for the 
lower boundary of the turbulent layer as: 

u(0,t) = akccos [k(x-ct) ] , 
w(0,t) = akcsin [k(x-ct) ] . 

In terms of the stream function the problem was solved 
with the boundary conditions: 

A(0,t) = ac, A (0,t) = -akc, (3.4) 

B ( 0 , t ) = B z (0,t) = 0 , (3.5) 

A (H , t ) = A z (H,t) = B (H , t ) = B z (H , t ) = 0 . (3.6) 
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Z (rn ) 




Figur e 1 . 



Undisturbed mean velocity profile (m/sec). 
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Consequently, the boundary layer was limited to a certain 
height H, where the wave-caused velocity perturbations were 
set equal to zero. 

Yefimov (1970) studied a similar problem in which the 
upper limit of the boundary layer was chosen sufficiently 
high so that it was reasonable to assume the conditions of 

(3.6) to hold. This upper limit was found to be greater 
than one and one-half wavelengths in order to reach a sta- 
tionary solution. 

C. RESULTANT FACTORS 

As a consequence of the initialization functions (3.1), 
(3.2) and (3.3); the boundary conditions (3.4), (3.5) and 

(3.6) , the functions A(z,t) and B(z,t) were time stepped 
forward until a steady state solution was reached. At each 
time step the Reynolds wave stress was computed from: 

T = -uw = - •jk (B^A - A z B) . (3.7) 

According to equation (2.17) a non-zero Reynolds stress will 

result in a mean horizontal velocity (u) of the wave-caused 

perturbations. This supplement to the initial undisturbed 

mean velocity may be computed as: 

H 

u = - i $ Tdz (3 * 8) 
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and applied as a correction to U at each time step by: 

U * z 

U(t) = — In (— + 1) + u . (3.9) 

rs 2 

o 

D. FINITE DIFFERENCING 

The finite differencing was accomplished through the 
use of centered finite difference schemes for all terms 
except those involving friction. These terms were computed 
at time (t-At) . Examples of the schemes are: 

A z = 2 Az ^ A K*1 ” A K-1^ ’ (3.10) 

and 

A zz = 2 ^ A K+1 - 2A K + A K-1^ ' (3.11) 

Az 

As a result of the finite difference schemes and the 
boundary conditions, the variables at the second grid point 
are specified by: 

B (Az , t ) = 0 , 

In addition the difference schemes and the boundary condi- 
tions specify that the stream function amplitudes A(z,t) 
and B(z,t) remain constant at both the bottom and top two 
grid points . 

E. NUMERICAL SOLUTION 

Equations (2,15) and (2.16) were solved utilizing the 
principle of Richtmyer (1957) with a forward time step to 
start the marching procedure. To avoid solution separation, 
the problem was restarted every 50 time steps with a finite 
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difference scheme developed by Matsuno (1966), This was a 



two-step iteration to simulate the Euler backward differ- 
ence method . 

All calculations were performed to double precision 
accuracy and continued until, the Reynolds stress (3.7) 
reached a steady state condition which is defined by: 



A space step of 0.25 meters and a time step of 0.0125 
seconds were utilized. Additional accuracy was attained 
in some cases by reducing the space step to 0.125 meters 
with a corresponding time step of 0.00312 seconds. 




o 



1x10 



-5 
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IV , COMPUTATIONAL PROCEDURES 



A. VARIABLES 

The forecast formulas, as stated in (2.15) and (2,16) , 
allow for five independent variables. These are wave 
amplitude (a), wave number (k) , dynamic velocity (U^) , rough- 
ness length (z ) | and the coefficient of turbulent exchange 
o 

(V) • 

Yefimov (1970) assumed that the slope of the waves was 
small, so that the expressions for the velocities on the 
boundary predicted by the theory of waves of small slope 
could be applied. This assumption of a small slope is 
applied by establishing the relationship ak — ,01, which is 

incorporated in the model as a means of computing the 
initial wave amplitude A(0,0), 

The resulting variables k, U^, Z Q , V were then con- 

sidered systematically by first determining the range of 
variation of each of the variables. Then central values 
were selected and used as example cases by holding three 
variables fixed and computing test cases of the remaining 
variable over its full range, 

Kraus (1967) studied the relationship of the roughness 
length to that of the mean wind where roughness length is 
given by : 
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The constant, m , was determined by different researchers 
with values ranging from 0.0125 to 0.08. A central value 
of m = 0.042 was selected and the roughness length was com- 
puted from: 

z = 0.00428 U 2 , (4.1) 

o * 

for the case where z was dependent upon U.. This reduces 

o x 

the independent variables to U^, k and V. 

The majority of the experiments were conducted with 

z « 0.0001 meters after Hill (1962). 
o 

B. MODEL CONSIDERATIONS 

Although the value of the coefficient of turbulent ex- 
change is unknown, it is reasonable to assume that it can 
be variable with height as discussed by Yefimov (1970). 

This possibility has been incorporated in the model and 
several cases have been investigated where V was allowed 
to vary with height. 
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V 



RESULTS 



A. EXPERIMENTAL RESULTS 

The number of variables and their large range of varia- 
tion necessitated numerous numerical experiments. To re- 
duce the number of experiments as much as possible, the 
variables were restricted to a limited range. The selection 
of a variation increment was guided by the necessity to 
choose a sufficiently large increment while remaining with- 
in the scope of full investigation of the problem. The 
variables and the range of values used are shown in Table 1. 
These were used in the various combinations as shown on 
Figures 2 - 10. 



TABLE I. 


Range of the variables 
numerical experiments . 


cons idered 


in the 


k (m 1 ) 


U A (m/ sec) 


V (m^/sec) 


z (m) 
0 


0.050 


0.010 


1.100 


0 . 10 


0.075 


*0.020 


*0 . 240 


0 . 01 


0.100 


0.040 


0.200 


0.001 


0.125 


0.060 


0.110 


0.0005 


0.150 


0.080 


0 .020 


*0 . 0001 


0.175 


*0 . 100 


0.011 


0 . 00005 


*0.200 


0.150 


0 . 002 


0 . 00001 


0:225 


*0 . 200 




0 .000005 


0.250 


0.250 




0.000001 


0.275 


0.300 






0.300 


0.350 

0.400 







^central values 
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The existence of a critical level was an important 
factor in the determination of the profile of the solution 
The critical level is defined as that height where phase 
speed is equal to the mean wind speed. Since the phase 
speed is inversely proportional to the wave number(c = 

/G 7k) , the phase speed will decrease as wave number in- 
creases. Therefore, for a fixed mean wind profile, an 
increase in wave number will decrease the height of the 
critical level (if existing) . In the cases of fixed wave 
numbers, an increase in U A , or decrease of Z q should de- 
crease the height of the critical level. 

The first case considered was that of a variable wave 
number. In this case the results were greatly affected 
by the value chosen for U A and therefore two experiments 
were conducted. In the first experiment (Figure 2) the 
wave number was varied with U*= 0.20 m/sec. The second 
experiment (Figure 3) was identical but U A = 0.020 m/sec. 

Some noteworthy features of these two tests are: 

a. the decrease in magnitude and variation of T as 
U A was decreased; 

b. the decrease in magnitude and variation of T as 
k was increased; 

c. a maximum value of T occurs at approximately one 
meter above the wave surface; 

d. the increase in T at a certain wave number when 
U A = 0.20 m/sec. It is significant that this was the 
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first case in which the critical level was within the layer 
where x varies significantly. 

The second case considered was that of different mean 
wind profiles. The mean wind given by (3.1) is dependent 
upon by U. and z . Since U. and z are considered to be 
interrelated there are numerous ways to vary the mean wind 
profile. In this study three experiments were conducted. 
The first case (Figure 4) was that of a fixed wave number, 
roughness length and turbulent coefficient. In the second 
case (Figure 3) the wave number and turbulent coefficient 
remained fixed but the roughness length was allowed to vary 
according to (4.1). Hence, both dynamic velocity and rough 
ness length varied in this case. The third case (Figure 6) 
was that of fixed dynamic velocity,' wave number and turbu- 
lent coefficient, while roughness length was allowed to 
vary. 



Some notable features are: 

a. the increase in magnitu 

stress as z. was decreased; 
o 

b. the increase in magnitu 
stress as was increased unti 
to 0.25 m/sec. This is the po 
is within the layer where T var 

The third case to be consid 
turbulent coefficients (Figure 
were held fixed. 



de of the wave-related 

de of the wave-related 
1 is approximately equal 

int where the critical level 
ies significantly, 
ered was that of various 
7) while the other variables 
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A noteworthy feature is the increase in magnitude of T 

2 

as V decreased to 0.020 m /sec. This is followed by a 
decrease in the magnitude of T as V continued to decrease. 

In the fourth case the turbulent coefficient was 
allowed to vary with height (Figure 8). In the first ex- 
periment the turbulent coefficient was decreased with 
height, to 10 meters, according to the relationship: 

V ( z ) = 0.5 sinh (H-z ) / s inh (H) . 

On the subsequent experiments V increased linearly with 

2 

height (Figure 8) until V = 0,24 m /sec and then remained 
constant thereafter. The results of the first experiment 
are shown in Figure 9 and those of the increasing V experi- 
ments are shown in Figure 10. A similar case, but with V 
held constant, is shown in Figure 10 for comparison. 

Some notable features are: 

a. the similarity in the profiles of the cases of 
increasing V with that of decreasing V; 

b. the presence of a maximum at approximately 0.5 
meters above the wave surface; 

c. the rapid change from positive to negative stress; 

d. the decrease in magnitude of T as V increases to 
greater heights. 

B, COMPARISONS WITH DATA 

Additional experiments were conducted to compare the 
results of this model with that of the data analyzed by 
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Davidson and Frank (1973) . The initial input parameters 

(Table II) were combined with a constant turbulent coeffi- 

2 

cient (V = 0.240 m /sec) and the calculations time stepped 
forward until a steady state was reached. The resulting 
stress profiles and amplitudes are shown in Figures 11-13. 



TABLE II. Dynamic velocity, wavelength and roughness 
length values determined by Davidson and 
Frank (1973) with corresponding wave numbers 
and critical levels. 



Case 


U A (m/sec) 


A(m) 


z (m) 
o 


k(m 1 ) 


critical 
level (m) 


1 . 


0.200 


35 


0 . 003280 


0.179 


43.50 


2. 


0.167 


23 


0 . 000079 


0.273 


29.00 


3. 


0.121 


23 


0.000130 


0.273 


none 



Appendix E outlines the procedure for obtaining data from 
the reported results of Davidson and Frank (1973) and also 
the procedures for interpretation of the numerical results. 
The elements for comparison have been displayed in 
Table III. 

A' significant feature of Table III is the high-degree 
of correlation of the magnitude of the velocities obtained 
from the variance spectra, phase amplitude and numerical 
amplitude procedures. An equally significant feature is 
the correlation of the magnitude of the stresses obtained 
by the velocity cospectra with those from the numerical 
calculations. It should be noted that the sign of the co- 
spectra stress is dependent upon the wave phase. The 
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differences in the signs should therefore not be considered 
as a significant discrepancy. 

The comparisons are the first of this kind for data 
from the natural regime. The agreement is significant. 
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TABLE III. Comparison of • the numerically derived results with 
statistically analyzed results from Davidson and 
Frank ( 1973 ). 
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Figure 2 . Distri bu t i o n of the wave — caused stress (l in 
m^/sec^) with height for = 0.2 m/sec, V — 

0.240 m^sec^, z = 0.0001m and k equal, respec- 
tively, to: ° 1) 0.050, 2) 0.075, 3) 0.100, 

4) 0.125, 5) 0.150, 6) 0.175, 7) 0.200, 

8) 0.225, 9) 0.250, 10) 0.275, 11) 0.300. The 
arrows indicate the height of the critical level 
and the identifying wave number. 
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Figure 3. Distribution of the wave-caused stress (T in 
m^sec^) with height for = 0.02 m/sec, V = 

0.240 m2/sec, z = 0.0001m and k equal, respec- 
tively, to: ° 1) 0.050, 2) 0.075, 3) 0.100, 

4) 0.125, 5) 0.150, 6) 0.175, 7) 0.200, 

8) 0.225, 9) 0.250, 10) 0.275, 11) 0.300. 
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Figure 4. Distribution of the wave-caused stress (t in 
m^/sec^) with height for k = 0 . 2 0 0m~ 1 , V = 

0.240 m^/sec, z = 0.0001m, and U* equal, respec- 
tively, to: ° 1) 0,01, 2) 0.02, 3) 0.04, 

4) 0.06, 5) 0.08, 6) 0.10, 7) 0.15, 8) 0.20, 

9) 0.25, 10) 0.30, 11) 0.35, 12) 0.40. The 
arrows indicate the height of the critical level 
and the identifying dynamic velocity. 
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Figure 5. Distribution of the wave-caused stress (t in 

m^/sec^) with height for k = 0.20m“l, V = 0.240 
m^/sec, z = 0.00428 and U^, equal, respective- 

ly, to: ° 1) 0.01, 2) 0.02, 3) 0.04, 4) 0.06, 

5) 0.08, 6) 0.10, 7) 0.15, 8) 0.20, 9) 0.25, 

10) 0.30, 11) 0.35, 12) 0.40. The arrows 
indicate the height of the critical level and 
the identifying dynamic velocity. 
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Figure 6. Distribution of the wave-caused stress (x in 
m^/sec^) with height for U = 0.100 m/sec, 
k = 0.200 m~l, V = 0.240m 2 ?sec and z equal, 
respectively, to: 1) 0.1, 2) 0.01, 3) 0.001, 

4) 0.0005, 5) 0.0001, 6) 0.00005, 7) 0.00001, 

8) 0.000005, 9) 0.000001. 
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Figure 7. Distribution of the wave-caused stress (t in 
m 2 /sec) with height for k = 0.200 m -1 , = 

0.100 m/sec, z = 0.0001m, and V equal, 
respectively, to: 1) 1.100, 2) 0.200, 

3) 0.110, 4) 0.020, 5) 0.011, 6) 0.002. 
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Figure 8. Variation of the turbulent coefficient with 
height, respectively, to: 1) decreasing, 

2) increasing to 5m, 3) increasing to 10m, 

4) increasing to 15m, 5) increasing to 20m, 

6) increasing to 25m. 
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Figure 9. Distribution of the wave-caused stress (x in 
m^/sec^) with height for = 0.100 m/sec, 

k - 0.200 m“l, z = 0.0001m and V decreasing 
with height. 
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Figure 10. Distribution of the wave-caused stress (t in 

m^/sec^Y with height for U^ = O.lOOm/sec, k = 

0.200m“^~, z = 0.0001m and V increasing with 

height, respectively, to: 1) 5m, 2) 10m, 3) 

15m, A) 20m, 5) 25m.- Curve 6) is the case of 

constant V with the same U., k and z values. 

* o 
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T >- 

figure 11. Distribution of the wave-caused stress (t in 
m^/sec^) with height for k = 0.179m“l, = 

0.200 m/sec, z = 0.003280m, and V = 0.240m z / 

2 ° 
sec . 
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Figure 12. Distribution of the wave-caused stress (t in 
m^/sec^) with height for k = 0.273m , U. - 

0.167m/sec, z = 0.000079m, and V = 0.240m / 

2 o 

sec'-. 
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Figure 13. 



Distribution of the wave-caused stress (l in 
m^/sec^) with height for k = 0.273m“-*, = 

0 . 121m/sec , z = 0.000130m, and v=0.240 m / 

,2 ° 



sec' 
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VI 



CONCLUSIONS 



The Reynolds wave stress of the perturbation-velocity 
field in the atmospheric near water layer above waves has 
been obtained from a numerical model. The resultant Rey- 
nolds stress can be used to compute the energy which is 
transferred by the average motion of the air to the wave 
fluctuations and vice versa. Thus a positive stress indi- 
cates the wind energy maintains the wave fluctuations sig- 
nifying a transfer of energy from the wind to the waves. 

The maximum Reynolds stress and greatest energy trans- 
fer was found to exist for the largest mean wind profile 

(largest U, and smallest z ), the lowest wave numbers and a 
* o 

2 

turbulent coefficient of approximately 0.020 m sec. In 
addition, the maximum stress occurred at approximately one 
meter above the wave surface. 

A striking result is the increase in the magnitude of 
the Reynolds stress as the mean wind velocity approaches 
that of the wave phase speed in the lowest levels of the 
boundary layer. For those cases where the critical level 
occurs in the region of the maximum stress fluctuations 
there is a significant increase in the magnitude of the 
stress. Thus the energy is greatest when wave velocity is 
close to the mean wind velocity in the lower levels. 

Another significant result is the change which occurred 
in the stress profile when the turbulent coefficient was 
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of this change 



allowed to vary with height. The net effect 
is a decrease in the total energy transfer. 

The comparison of numerically generated velocities and 
stresses with those of measured data yielded a very high 
degree of correlation. These results will justify future 
numerical research and data collection to provide new in- 
sights into the air-sea energy transfers. Specifically, the 
correlation of a numerical model with measured wind data 
could determine the characteristics of the turbulent coeffi- 
cient due to the wave-caused motion. Concurrent studies 
being conducted by Aurand (1973) and Strieker (1973) may be 
considered as sources for additional data to correlate the 
statistically and numerically derived stress patterns. 
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APPENDIX A 



BAROTROPIC INSTABILITY TEST 

Neuman (1969) considered a simple barotropic numerical 
model. The object of the test v as to determine the growth 
rate as a function of the wave number. The similarity be- 
tween the present study and that of Neuman allowed for a 
check of the present model for the special case of no 
f r ic t ion . 

However, to accomplish the check, several modifications 
were made to the present model. These are: 

A(z,0) = cos (“^-) , and 

n. 

V ( z ) = 0 . 

The mean wind velocity profile is given by: 

U ( z ) = -U tanh (y/y ) (Al) 

o o 

as shown in figure 14. 

The procedure required marching the calculations for- 
ward in time until the growth rate stabilized. The point . 
of maximum Reynolds stress for three different times were 
then compared. The magnitude of the stress at these points 
was found to be growing exponentially, as was expected. 

As an additional check, the growth rate of the Reynolds 
stress was evaluated and compared to the growth rate as 
given by the function: 

2U n 

exp( — ) t . (A2) 

y o 



45 



Z(m) 




Figure 14. Mean wind profile for U - -U q tanh(y/y Q ), 
(m/ sec). 



4 6 



In this case the stress at the maximum point was equated 
to (A2) and the equation solved for n using: 

U =9.0 m/sec , 
o 

y = 4 . 0m . 

The resultant growth rate was 0.181. This is more 
accurate than that of Neuman who found a growth rate of 
0.160, but it is less than that of Betchov and Criminale 
(1967) who set the rate at n = 0.19, The latter difference 
is because the present solution has been solved with a 
finite boundary. 
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APPENDIX B 



EXAMINATION OF THE VISCOSITY TERMS 

In this test the initial field of the cosine stream 
function amplitudes was computed utilizing equation (3.2) 
At this point the mean wind and phase speed were set to 
zero. Then the forecast equations were marched forward 
until a steady state was reached. 

These numerical values were then compared with the 
analytical values derived from the solution of the equa- 
tion : 



A ■ + 2A + A = 0 . (Bl) 

z z z z z z 

The boundary conditions utilized in the solution of (Bl) 
are similar to those given by (3.4) and (3.6). These 
conditions are: 

A(0) = 2c , A z (0) = -2kc , 

A (H) = A (H) = 0 . 
z 

The solution of (Bl) is given by: 

A(z) = c^ sin(kz) + c^ cos(kz) + c^ kz sin(kz) 

+ c^ kz cos(kz) , (B2) 
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where : 



C 1 " r IA Z (0) - c 4 k] ’ 

c 2 = A , 

C 3 = ' Hk si'n"(Hk ) t A z ( 0) s in (Hk) + A(0)cos(Hk) + 



sin(Hk)+ Hk cos (Hk) ] , 



c , = 



H sin (Hk) 
:Hk)- 
• A (0) 



sin^ (Hk) 



hi l 

2 2 V 

-H k > 



[A(0)+ 



A (0) 



Hk 



-] s in (Hk) 



+ I- 



H 



+ A(0)k cot (Hk) ] cos(Hk) ] 



)■ 



The comparison of the numerical results with the 
analytical results derived from (B2) indicated that the 
minimum occurred at the same height for both cases. How- 
ever, the magnitude of the numerical minimum was smaller 
than that of the analytical results. 
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APPENDIX C 



EXAMINATION OF THE CASE OF NO MEAN WIND 

Yefimov and Pososhkov (1970) derived an analytic solu 
tion to the forecast equations (2.15) and (2.16) for the 
special case U = 0. This solution is given by: 



-u 1 w 1 = 



-2zk . ~zk//Za ,zk N x 

-ae + ae cos ( ) , (Cl) 



/2a 



where : 




c 



The condition U = 0 was established and the numerical 
solutions of (2.15) and (2.16) were computed to a steady 
state condition. The results of the analytical and numeri 
cal calculations of stress were compared for several dif- 
ferent wave numbers and for different values of the turbu- 
lent exchange coefficient. 

The results of. the stress obtained by both methods 
were very similar. In all cases compared the analytical 
stress reached a minimum at a lower height and somewhat 
smaller magnitude than the numerically derived stress. 

This discrepancy, however, has been primarily attributed 
to the scale of the grid increment. 
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APPENDIX D 



EXAMINATION OF THE RICHTMYER SCHEME 



This test was devised solely to test the method of 

solving equations as given by Richtmyer (1957). The test 

was conducted by equating the right hand side terms of 

TT Z 

equation (2.15) to sin ( ) where w = (H-3)Az. The result 

w 

is : 



dz- 



, 2. 3A . 

- k ) ^ - sin <— ) , 



and upon solution yields: 



A 



-1 

TT 2 2 

[(“) + k ] 

w 



At . 



(Dl) 



The numerical solu 


time step. 


Then the m 


amplitude , 


A(z , t) , was 


from (Dl) . 


The result 


to within 


three per cen 



tion was marched forward one complete 
inimum value of the stream function 
compared with the value derived 
was a comparison of the values 
t of the difference. 
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APPENDIX E 



EVALUATION OF ANALYZED DATA 

Davidson, and Frank (1973) reported observed values of 
wave heights and velocity fluctuations at two levels above 
the waves. The results from their spectral and joint pro- 
bability density function, conditional mean function (JPDF 
CMF) analyses are shown in Figures 15 and 16. 

The velocity variance spectra graphs show plots of 
n • <j) (n) vs. ln(n) where n is the frequency of the wave in 
Hz. The area under the significant spikes of the curve 
(shaded areas of Figure 15) can be used to determine the 

value of the square of the wave-caused velocity fluctua- 
2 2 

tion (u 1 or w T ). These areas may be computed by: 

Area = n • <£(n) * A ln(n) , (El) 

where A ln(n) is defined at the significant peak. 

The velocity cospectra graphs (Figure 15) are similar 
to those above but the area in this case is indicative of 
the wave-caused Reynolds stress (u'w 1 ). These results 
have been presented in Table III for comparison. 

Davidson and Frank (1973) also presented phase-ampli- 
tude results for the same data as above. These results 
(Figure 16), which may be interpreted directly, are shown 
in Table III for comparison. 
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The numerical maximum amplitude has been computed 
from the steady state values of the stream function ampli- 



tude by: 




(E 2) 



and 



, 1 / a 2 J . - a 2 \ 1 / 2 

w = k(A + B ) 



(E3) 



which follow as a result from equations (2.11) and (2.14). 
The numerical results were computed with the small wave 
slope assumption where ak — 0.1. However, this exact con- 
dition is not normally met in nature. Therefore, the 
values obtained from equations (E2) and (E3) , as well as 
the numerically derived Reynolds stress, must be scaled 
to reflect the actual wave slope. 

The calculation of a scaling factor first involved the 
determination of the actual amplitude (a^) . This may be 
interpreted directly from the "waves" diagram in Figure 16. 
The amplitude used for the numerical calculations may be 
computed by: 



The resulting velocity scaling factor (Fs) may now be ex- 
pressed by: 



a 



D 



Fs 



a 



9 



(E4) 



N 



and the quadratic stress scaling factor (Fq) by: 



2 



a 



Fq 



2 



(E5 ) 



a 



N 
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VELOCITY VARIANCE SPECTRA n.*W 
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LARGE WAVES . , LARGE WAVES / LARGE WAVES 
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Figure 16. Phase-amplitude results for large waves taken from. Davidson and 
Fr ank (19 7 3) . 
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